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Four-quark correlators and the factorization hypothesis are analyzed in the meson sector within 
Chiral Perturbation Theory. We define the four-quark condensate as Vvaix^o {T{qq) {x) {qq) (0)) , which 
is equivalent to other definitions commonly used in the literature. Factorization of the four-quark 
condensate holds to leading and next to leading order. However at next to next to leading order, a 
term with a nontrivial space-time dependence in the four-quark correlator yields a divergent four- 
quark condensate, whereas the two-quark condensate and the scalar susceptibility are finite. Such 
a non-factorization term vanishes only in the chiral limit. We also comment on how factorization 
still holds in the large Nc limit, provided such a limit is taken before renormalization. 
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I. INTRODUCTION 

Scalar condensates play a relevant role in QCD, since they are directly related to vacuum properties. The quark 

condensate (qq) is a parameter deeply related to spontaneous chiral symmetry breaking and the description of low- 

-^ ^ energy QCD. In principle, quark condensates of arbitrary order {{qq)"') are also built out of chiral noninvariant 

II ■ operators with the vacuum quantum numbers and are also related to chiral symmetry restoration. In addition, quark 

CIh' condensates appear directly in QCD sum rules, through the operator product expansion (OPE) approach 1], vifhere 

^ . the following hypothesis of factorization or vacuum saturation is customarily made: 

>. 

^^ , Note that we have particularized to the case where the four-quark operator has the quantum numbers of the scalar, 

[ ' isoscalar, colorless condensates that we are interested in. In addition, N = iNcNj, where Nc and Nf denote the 

^_J ' number of colors and flavors respectively and g is a Dirac spinor, flavor and color vector. We remark that in the 

. • large- Nc limit factorization simply reduces to {{qq)'^) = {qq)"^- The second term in Eq.([l]) comes from the contraction 

'^ ' of indices (including color) between the first and second qq operators. 
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The use of the factorization hypothesis is a key point in order to estimate the size of higher order condensates in the 



^^ ' OPE. However, its justification is still a matter of debate. It was shown in y that factorization implies that {{qq)^) 
becomes dependent on the QCD renormalization scale. This means that for QCD sum rules including six-dimensional 
operators, like (qq)'^, one cannot write a renormalization-group (RG) invariant four-quark condensate, preventing RG 
improvements of such sum rules. This is not a problem when considering six-dimensional pure-gluon operators or 
quark operators with dimensions lower than six, like the RG-invariant qA4q with A4 the mass matrix. We will come 
Cd ' back to this point in section HVl The validity of vacuum saturation has also been questioned within the framework of 
finite-energy sum- rules [3] and has been formally shown not to hold when dressed QCD vertices are considered Q. 

In this work we will present a study of the scalar four-quark condensate within the framework of Chiral Perturbation 
Theory (ChPT). Since ChPT relies only on symmetries and not on vacuum saturation or dominance assumptions, as in 
some of the approaches commented above, it will allow us to obtain low-energy model-independent results concerning 
the factorization hypothesis. 

An important point concerns the definition of the quark condensate in terms of Green functions. In the chiral 
lagrangian framework, one has not access to individual quark operators at a given space-time point a;, but to the 
low-energy representation of the quark-antiquark operator qq{x), given by a functional derivative with respect to an 
external scalar source (see details in section |TT1) . Therefore, a natural way to define the four-quark condensate is 
through the limit of the two-point function (four-quark correlator) : 
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{{qq)') = ]im{T{qq)ix){qqm) (2) 

This is the definition that we will choose to work with here, where all the divergencies will be treated within the MS 
scheme in dimensional regularization, as it is customary in ChPT. However, from the comments above, it is not clear 
that the four-quark condensate itself has to be a scale-independent and finite object, which means that the a; — s- 
limit is ill-defined and other definitions in terms of Green functions could give different answers. Actually, Eq.([2]) is 
not the usual MS definition when working for instance with four-quark vacuum expectation values in the context of 
electroweak penguin contributions [5|, |6|, where the following prescription is used instead: 

{{qqf)= j d^x{T{qq){x){qqm)5^^Hx)= j ^n{Q^) (3) 

where the integrals are defined in Euclidean space-time dimension D and I{[Q^) is the Fourier transform of the 
correlator {T{qq){x){qq){0)). In the ChPT framework, we will show (details are given in Appendix [Bt that this 
definition gives the same result as the one in Eq.Q meaning that factorization is spoiled at next to next to leading 
order (NNLO), which questions seriously the validity of the factorization hypothesis, now from the point of view of 
the low-energy representation. 

The four-quark two-point correlator, apart from defining the four-quark condensate, is also related the chiral or 
scalar susceptibility, defined as x = -~d{qq)/dmq and which can be written also in terms of {T{qq){x){qq){0)). The 
susceptibility is a crucial observable regarding chiral symmetry restoration, since it is associated to thermal fluctuations 
and tends to grow near the critical point p]. For us, the susceptibility will serve as a crucial consistency check, since 
we can calculate it directly as a quark mass derivative or through the four-quark correlator and both should coincide 
and be finite and scale independent. 

Therefore, we will give the complete results in ChPT for the four-quark correlators and four-quark condensates in 
SU{2) and SU{3) up to NNLO, performing a consistency check by calculating the scalar susceptibility and showing 
the robustness of the result under different definitions of the vacuum four-quark expectation value. In addition, 
the discussion of factorization breaking necessarily implies the calculation and renormalization of the two-quark 
condensate also at NNLO, which we will perform also explicitly here. We will also carry out the large- iVc analysis of 
the factorization breaking, which can also be performed from the low-energy representation and is formally relevant. 
These are the main results of this work. 

The plan of the paper is the following: In section |lT] we present our calculation of the relevant four-quark correlators 
for two and three flavors. The details of the calculation are given for Nf = 2 for simplicity. The scalar susceptibility 
derived from the four-quark condensate is obtained in section IIIII The factorization hypothesis is then examined 
in section IIVI whereas in section |V] we discuss the large- TVc limit of our results, regarding factorization. In section 
IVIl we present a brief summary and our conclusions. Finally, in Appendix \X\ we provide the detailed mathematical 
expressions for the two-quark condensates to NNLO in ChPT and discuss in detail their renormalization, whereas 
in Appendix [B] we show the equivalence of our deflnition of the four-quark condensate with the usual one in the 
literature. 

II. FOUR-QUARK CORRELATORS 

Our main object of study wiU be the time-ordered four-quark correlator {T{qq){x){qq){0)). We will follow the 
external source method and write this four-quark correlator as a second functional derivative of the QCD generating 
functional Zqcd[s] with respect to the scalar source s{x), which in general will be a matrix- valued function in flavor 
space and couples to the QCD Lagrangian as: 

Zqcd[s] = 'Dq'Dq...expi d'^xjCQCD[q,q,s{x),...], 

^QCd[s] = q{ip~s{x))q+..., (4) 

where the rest of the Lagrangian and fields indicated by dots are irrelevant for our purposes and sum over Nf light 
flavors, Nc colors and Dirac indices is assumed in qq. The physical QCD Lagrangian and partition function correspond 
to setting s{x) — A4, the quark mass matrix, in the above equation. 



We will consider the effective low-energy representation of Zqcd [s] given by Chiral Perturbation Theory [8| , built 
from chiral symmetry invariance as an expansion in external momenta (derivatives) and meson masses: 

Zqcd[s] ^ Zeff[s] = / Vcfy^expi / d'^xCeff[4>'^,s{x)], 

^eff ^ C2 + Ci + Ce . . . , (5) 

where the subscript in the effective Lagrangian indicate the order in the derivative and mass expansion, formally 
Ck ~ 0{p^) [s — 0{p'^) in the standard ChPT power counting). Note that (j)"- denote the NGB fields, usually 
collected in the SU{Nf) matrix U = exp[«Aa(/>"/i^], where Xa are the Gell-Mann or Pauli matrices for Nf = 3 and 
Nf = 2, respectively and F is the pion decay constant in the chiral limit. The Lagrangian £2 is the non-linear sigma 
model: 

/:2 = ^Tr[9^f/t9^C/ + x(C/ + C/t)], (6) 

with X = 2Bos{x). When s{x) — Ad, the constants mq,F,Bo appearing in £2 are related to meson masses, decay 
constants and to the quark condensate. For simplicity, we will work in the isospin limit ruu = rud = m, so that, to 
lowest order in SU{2), M§^ = 2mBo{l + 0{p^)), F^ = F(l + 0{p^)) and (qq) = BqF{1 + 0{p^)). As usual Mo^,oif,o^ 
stand for the leading order meson masses, in terms of which we will express our results. Their relation to the physical 
masses is given in Eas. (|A9p and (lAlOp in the Appendix \K\ In addition, and for our purposes here, Weinberg's chiral 
power counting 9], on which Chiral Perturbation Theory relies, can be equivalently accounted for by keeping trace 
of inverse powers of F , which will be used extensively along this work. 

The Lagrangians £4 and Cq are given in Q and [ifll, respectively, where use has been made of different operator 
identities, partial integration and the equations of motion to the relevant order. Those lagrangians contain the so- 
called low-energy constants (LEG) multiplying each of the independent terms compatible with the symmetries. The 
£4 LEG receive different names depending on whether they multiply terms containing U fields or not; respectively, 
Li and Hi in the SU(3) case. The terms without U fields are contact terms containing just external sources and no 
fields, but they are needed to absorb some divergences coming from loop diagrams using £2 vertices. The original 
SU(2) lagrangians in [ll| are written in terms of vector fields instead of matrix fields U as above, but they also use 
different names for the £4 low-energy constants - U and hi in this case. However, it is possible to recast [IJ] these 
lagrangians using matrix field notation, that we will use throughout this paper, and keep the same L, /li low energy 
constants. The relation between the SU(3) and SU(2) low-energy constants is given in [8|, [i3| and |14l |. 

This name differentiation for the C% is not followed any longer [lOj: all of them are called Ci in the SU(2) case 
and d in the SU(3) case. Note that the 0{p^) LEG contained in £5 absorb both two-loop divergences from £2 and 
one- loop divergences in diagrams with £4 vertices. All the details for renormalization of quark condensates up to the 
order we are considering here are given in Appendix \K\ We recall that the £4 Lagrangian in SU{2>) contains also 
the Wess-Zumino-Witten (WZW) [la| anomalous term, accounting for anomalous NGB processes, whose coefficient 
is fixed by topology arguments and is proportional to the number of colors Nc- 

A. Two flavors 

For simplicity, we will discuss the full details of our approach in the simpler case Nf — 2. Thus we will denote by 
the subscript I the light quark correlator, and study {qq)i = uu + dd . Note that we have defined the scalar source 
s{x) as a matrix, but since for the physical partition function it corresponds to the mass matrix A4, which is diagonal, 
we are thus only interested in the diagonal elements of s{x) and we can set the rest of the source terms to zero. In 
particular, for the two flavor case Ai — ml2 and we can write s{x) = sq{x)12, so that: 
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Proceeding in the same way now for the four hght quark correlator: 

-1 A A 

(T{qq)lix){qq)im = -^^^ , , ., m^ ^-ff^'^^ 

Zeff[m\ 6so{x) dso(O) 



5sq{xY I^ \ 8sa(x) (5so(0) 



We will regularize all our expressions in dimensional regularization with Z? = 4 — e and for that purpose we keep the 
i'-dependence in the (5-function term above. 

Now, from Eq.([8]), and using the Lagrangians in [3,[lO], we obtain the following result: 
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{T{qq)i{x){qq)i{0))NNLO = {T{qq)i{x){qq)imNLO + ^B^F^ 
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where the NNLO constants Ci arc defined in Ea. (|A3p and, as usual 
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Note we have defined K^^\x) as the connected part of the four-pion correlator to leading order: 



K^^\x) = {Trix)M^)^\0)M0))LO - (T0"(O)0,(O))io = QGlix), 
G.(a:) being the pion propagator to leading order and the factor of 6 = 2(A^? 



)LO ~ [-1-9 y^)(Payy}))LO = ou^i^a;;, (12) 

J'i — 1) comes from the Wick contractions 
and is nothing but twice the number of NGB fields. The details of the renormalization and the dependence of the 
constants II, h\ and c[ on the renormalization scale jjl are given in Appendix \K\ 

To understand the structure of the different contributions to Eqs.® and ((TU)) it is useful to recall the general form 
of the SU{2) low-energy Lagrangian terms depending on the external scalar source. For our NNLO calculation, we 
will need to keep terms up to 0{F^^). Let us then separate the terms in the Lagrangian [a, [i3|) according to their 
s-dependence after expanding the U in NGB fields: 
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(13) 



where we have also made explicit the leading l/F^ dependence of each term. The superscripts "nt/)" indicate the 
number of NGB fields or field derivatives on each Lagrangian contribution. Note that, since £fe = 0(j>^) in derivatives 
or s-powers (s = ©(p^)), it counts at least as ©(l/F*^"^), but the XjF"^ order of each term grows when increasing 
the number of NGB fields, (/>. We have represented the vertices arising from the different pieces of the Lagrangian 
above on the left column of FigHJ Note that all £5 terms in Eg. p^ have the 00 superscript, because, to this order, 
they are simply constants. The constant L^ term enters in {(l^tlf n j^ lq ^^'^ ensures that one can renormalize the full 
result so that the quark condensate is finite and scale- independent. The term containing (9so)^ does not contribute 
to this order. The details as well as the explicit expression of the condensates up to NNLO are given in Appendix \K\ 

Once the structure of the vertices arising from the Lagrangian Ea.(|13p are understood, we represent diagram- 
matically in Figure [T] the different contributions to {T{qq){x){qq){0)). On each diagram, the horizontal dotted line 
represents spacetime, where each quark antiquark-bilinear stands at separate points and x. To LO and NLO - 
respectively 0{F'^) and 0{F^) - all contributions are disconnected, as seen in diagrams a, b and c. The reason is that 
we can only use the £3 term once and therefore, the NGB line has to close upon itself- a tadpole. This gives diagram 
(b) in Figure [1] To NNLO (©(F")) we have all the possibilities shown in Figure [1] in diagrams (d)-(j). If one of the 
vertices comes from £4 or Cq, once more there is at most one NGB line and the resulting diagram is disconnected. 
Note that among these is the ^'■■^•'(x) term in Eg. dTOl) from diagram (h). With only £2 vertices, one has a diagram 
with a double tadpole in one of the vertices leading to a LO propagator squared at the same point (diagram (d)), 
two vertices with one tadpole each (diagram (e)), a diagram like (b) but with the propagator renormalized to NLO 
(diagram (f)) and another with two NGB lines on each vertex but joined to form a connected one- loop diagram, which 
is diagram (j). Actually, the latter is the only possible connected contribution to this order, and gives the G'^{x) term 
in Eq. p^ . This whole discussion of vertices and diagrams will be valid also for the SU{3) case discussed below 

Let us now turn to the factorization hypothesis and the relation between the four quark correlation function and 
the two-quark condensate. We have collected in Appendix [X] all the two-quark condensates ChPT expressions up to 
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FIG. 1: On the left column we provide the diagrammatic representation of the vertices coming from the different terms of the 
Lagrangian in EQ. p3|l . The numbers attached to every vertex indicate the order of the Lagrangian. Diagrams a) to j) represent 
the different contributions to the four-quark correlator. The dotted horizontal line represents the spacetime separation between 
and X. Note that each NGB line decreases the order of the diagram by 1/F^ . Diagram j) is the first factorization-breaking 
term. 

NLO (given also in [Tl| for SU{2) and in ^ for SU{3)) and up to NNLO, which have been given explicitly in [1^] for 
SU{3). Numerical estimations including NNLO corrections are given in [ia)[l3- In view of Eqs.([7|).( jAlip and (jA12p . 
it easy to check that 



{T{qq)i{x){qq)i{0))NLO = {{qq}f)j^^Q, 

{T{qq)i{x){qq)imNNLO = {{qq)f) ^^.^^ + B^ Ustih + h,)6^''\x) + K^'\x) 



(14) 



We see that all contributions from disconnected diagrams in Figure [TJ other than the 5^-^^ term, can be absorbed in 
the two-quark condensate. Actually, up to NLO, we observe that {T{qq)i{x){qq)i{0)) in Eq.® is constant and equal 
to the NLO of the quark condensate squared, which leads to factorization in the iVc -^> oo limit (see section ITVl and 
W} . However, to NNLO the previous expression for a; = contains the G^(0) divergent contribution, even after the 
quark condensate has been renormalized and the 6^^' term regularized. We will show below that divergences cancel in 
physical quantities such as the scalar susceptibility, which is directly expressed in terms of observable quantities such 
as the free energy density. That is not the case for the four-quark condensate, which will remain divergent. Before 
analyzing these issues, let us extend the previous analysis to the SU{3) case. 



B. Three flavors 



In the SU{3) case, qq = uu + dd + ss, A4 = diag(TO, m, rus), s{x) — diag[so(a;), so{x), Ss{x)] and: 
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The s-dependent terms in the SU{3) effective Lagrangian are now the generahzation of Ea. (ll3p to include Ss{x), so 
that we have crossed terms hke SqSs, SqSs and so on, but the general structure is the same. As in the SU{2) case, 
the derivative terms (ds)"^ do not contribute to {T{qq){x){qq){0)) and thus only four Cq constant terms contribute to 
renormalization. As seen in Appendix \X[ they are proportional to the C'i LEG given in Eq. (jA3[) . Since we already 
presented the detailed discussion for the SU{2) case in the previous section, for the sake of brevity we cast our SU{3) 
results for {T{qq){x){qq){0)), which are much longer than before, directly in terms of the two-quark condensates, 
namely, 



{T{qq){x){qqm)NLO = {{qqf 
{Tiqq)ix)iqq)iO))NNLO = {{qq? 
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where K{x) is the extension of Eq. (|T2|) to the SU{2>) case: 



K{x) = {Tr{x)d^a{x)4>\^)U^))LO - {Tr{0)MO))lo = 2 [^GHx) + 4Glix) + Gl{x)] 



(17) 



(18) 



The ChPT expressions for the four-quark condensates to NNLO given in Eqs. p^ and pT|) (simplified in terms of the 
explicit expressions for {qq)NNLOi which are given in Appendix E)) . are among the main results of the present work. 
Note that, as it happened in the SU{2) case, the contribution Eg. pSj) stems from 2(iV? — 1) NGB propagators, 
although this time they have different masses. Similarly, we can calculate separately the strange and non-strange 
four-quark condensates, which also factorize up to NLO, whereas to NNLO we get: 
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(20) 



Once again the explicit expressions for the renormalized {qq)NNLO ar^ give n in Appendix Rl 

The {T{qq){x){ss){0)) correlator has been calculated up to NNLO in [I^, [3 i^ terms of the basis of the solutions 

to the Muskhelishvili-Omnes equations . 

We remark that the four-quark correlators to NNLO given in Egs. dT^ and (fTT)) are key ingredients to define the 

four-quark condensate and study the factorization hypothesis, as explained in the introduction. 



III. THE SCALAR SUSCEPTIBILITY 



In this section we will provide a consistency check of our calculation by analyzing the chiral or scalar susceptibility 
to the first nontrivial order, which can be obtained either differentiating the two-quark condensates or by integration 
of the four-quark ones. The susceptibility is defined in Euclidean space-time as 



XI = --^—{qq)i 

dm 



(21) 



and measures the condensate thermal fiuctuations, growing dramatically near the chiral restoration, as confirmed by 
different lattice studies 01. Therefore, let us consider the Euclidean (imaginary-time t = —it) version of Eq.(|3|) and 

-, — ) metric in the Lagrangian. Recall that the finite 
would correspond to r G [0, /3] with /? — 1/T. In addition, in 



dS]), replacing ijd'^x — >■ J dr J d'^x = J^d^x and the 
temperature T case, which we will analyze elsewhere 
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Eq.® and (fT6|) we have to replace —iS^{x) — )■ i5(r)(5^^ ^H^) = ^si^)- With these replacements, we can now relate 
the susceptibility with the four-quark correlators in the non-strange sector: 
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where Ve — /g d^x is the Z3-dimensional Euclidean volume and Z = Z[s = A4] ~ e"^^^ is the partition function 
with z the free energy density. 

The relation in Eq. (j22p between xi ^nd the four-point function allows us to check our previous results. From 
Eqs. (fT4l) and (fTQ]) . taking into account that: 



d^x[G.ix)]'^~^^Gm, 



(23) 



and the expressions Eqs. (IAll )- (jA2l) . together with the renormalization of the LEG in Eqs. dASI) and (jA6l) . we obtain 
using the last integral in Eg. ([2^ : 
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with Vi given in Eq.(| 

This is the same result that we get by taking directly the mass derivative of the quark condensate to NLO in 
Eq. (lAlip and (JAISP using the leading order relations between meson and quark masses [8]] . This represents a check of 
consistency of our calculation of the four-quark condensates to NNLO. In addition, we have explicitly checked (using 
Eq. (|A5p ) that the susceptibilities above are finite and independent of the scale /i. Furthermore, with the conversion 
between the SU{2) and SU{3) LEG given in [8]: 



lli^i) + hli^i) = 2 (^8Llif,)+2Lli^i)+H^,i^l) - -^uk - ^^„) 
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we end up with: 
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which is also consistent since the SU{'i) susceptibility is given by constant plus logarithmic terms in the mg -^ 00 
expansion, with no subleading terms in that expansion and therefore the very same expression has to be exactly 
recovered by calculating directly in the SU{2) limit. Note also that the susceptibility to this order is independent of 
F. 

Our result for the susceptibility is also consistent with a previous work [2l| where only the leading infrared order 
in the chiral limit was calculated, namely the logMg^ term inside the i/,r in Eq. (l24l) . This is the expected behaviour 
of the susceptibility from the 0(4) model universality class near the chiral limit and below the critical temperature, 
namely x ~ login, with mi the mass of the non-strange quark [7|, 121] . 

We can follow the same procedure to obtain the strange quark susceptibility in terms of our strange four-quark 
correlation function: 
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which, from Eq. ((20|) gives 
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We have explicitly double checked this result by taking the derivative with respect to mg of the NLO strange quark 
condensate in Eq. dAlSp . We remark that the results in Eqs. d^^ . (|25p and (l^5|) for the GhPT scalar susceptibilities 
have not been given elsewhere. 



IV. NON-FACTORIZATION 

As explained in the introduction, we define the four quark condensate through Eq.Q, although in Appendix IB] we 
show that this is equivalent to the more usual definition of Eq.®. Therefore, by taking the x — ;■ limit in Egs. p^ 
and (|17|) . and despite 5^^'{0) vanishes identically in dimensional regularization [22i] (now we are not integrating over 
X as for the scalar susceptibility), there is still a term that clearly breaks factorization, as defined in Eq.(IT]). In 
particular, we get in SU{2), from Ea. dTil) : 



(gg)2 ' F4 '^^ ' ' \F^ 

whereas in SU{3) from Eq.([T8]), we find: 



l + ^Gii0) + O -rl (29) 



^^'^'^^'^ -' ' 4l [3G'^(0) + 4G?,(0) + G^(0)] +o(^), (30) 



where the propagators G'i(O) are given in dimensional regularization in Eg. dAll) . 

The non- factorization terms above are divergent and independent of the LEG, once the two-quark condensate (qq) 
has been rendered finite with the renormalization of the 0{p'^) and 0{p^) LEG (see Appendix E^ . The renormaliz- 
ability of {qq) is of course consistent with the fact that qM-q is a QGD genuine observable RG-invariant. Therefore, 
our non- factorization GhPT results in Eq. ((29|) and (l30l) imply that the four-condensate is divergent and therefore the 
vacuum expectation value of (qq)'^ does not admit a meaningful low-j^nergy representation. 

Our result is consistent with the one-loop QGD RG analysis in [^, where only one flavor is considered. In that 
paper it is shown that factorization is incompatible with the renormalization group. Their argument goes as follows: 
the operator (qq)'^ mixes under renormalization with other four-quark operators, which can be chosen in combinations 
such that their vacuum expectation values would vanish if factorization holds. Then, assuming factorization for those 
other operators leads to the conclusion that {{qq)'^) is divergent, which in particular means that it does not factorize 
in terms of {qq)'^ and that one cannot write any RG invariant made of four-quark operators. 

Another interesting comment is that the factorization breaking terms in Eqs. (|29p and (|30p vanish exactly in the 
chiral limit, since then all dimensionally regularized propagators G7r(0) = Gk{0) — G^(0) = 0. In that case, we 
would be forced to examine the neglected NNNLO contributions in order to check the validity of factorization and the 
finiteness of the four-quark condensate. Recall that the arguments in j2] regarding four-quark operators hold actually 
for m = 0. 

V. LARGE Nc 

Let us now discuss the Nf and Nc dependence for the regularized expression, namely, before taking the D — 4 
limit. As we have checked for the SU{2) and SU{3) case, the connected four-field functions in Egs. p^ and ^TE\\ 
K{x) = 0{Ngb) = 0{N'i), where Nqb = Nj — 1 is the number of Goldstone Bosons. In addition, the Nc leading 
behavior of the different GhPT constants is well known Q from the QGD 1/Nc expansion. In particular, F^ — 0{Nc). 
Therefore, the first term that breaks factorization in Eas. (l29)) and (jSOj is O ( N'i /N^ ) , which is rather different from 

the \l{ANfNc) scaling suggested in Eq.([T]). Unfortunately, we cannot say much more about the Nf behavior of higher 
order terms, that could change the global Nf behavior. Note that the Nf dependence of the quark correlators has 
been studied in detail in [23i] with a different motivation. 

In the following, we will easily deduce the 1/Nc behavior and, in particular, we can study the large Nc limit before 
renormalization, and we will see that, in such formal case factorization holds for Nc — > oo. First of all, contrary to 
Eq.(IT]), in Eqs. d^Hl) and ([50]) there are no 0{\/Nc) terms. These could have arisen from contributions of the type 
LiG(fi) / F^ , when Li is 0{Nc), that actually appear in the calculation. However, as we have said before, the whole 
Li dependence of the four-quark condensate is exactly that of the two-quark condensate squared and thus such terms 
do not break factorization. The same happens with the 0{p^) Ci LEG in Ea. (jA3l) . Still, one could wonder if 0{1/Nc) 
or larger Nc powers could arise from higher chiral orders that we have not calculated explicitly here. 

Of course, as seen in Egs. iP^ and ([501) . these higher chiral orders count at least as 0{1/F^). Since F^ — 0{Nc), 
this already introduces a l/Nc factor, but it is not the only one, since the LEG can carry their own Nc behavior. In 
particular, we recall that, according to the chiral power counting discussed in scctlTTl the ©(l/F") contribution to 
the ratios in Egs. iP^ and (I50|) comes from connected diagrams with n = 2{L + 1) + J2d ^d{d — 2), with L the number 
of loops and Nd the number of vertices from Cd — C,2-,C^^.... Note that a non- factorizing term requires at least £ = 1, 



the leading contribution being the connected one-loop diagram (j) in Figure [T] with two £2 vertices. This diagram 
yields the factorization breaking terms in Egs. iP^ and ([50)1 . 

Now, the highest Nc scaling of the LEG from Cd is 0{Nc ). The reason is that these LEG, when divided by 

F''"'* should yield 0{Nc) contributions at most, as expected from the large- A'c behaviour of the low-energy generating 
functional [Sj. This includes the WZW term, which is the anomalous part of £4 and is multiplied explicitly by Nc 
[l5{ . Although the WZW term does not depend on the quark mass, it could enter in this calculation through loop 
contributions. It is possible, of course, that some LEG do scale with a smaller Nc power. For instance, the Li to 
LiQ appearing in £4, are known to scale as 0{Nc), except L4, Lq and L7, which scale as 0(1). These are model- 
independent QGD predictions obtained in Q, with the exception of iy, which was taken there as 0{Nc). This Ly 
counting corresponds to integrating the rj' as a heavy particle but then considering mi, ~ 0{1/Nc) and therefore a 
light particle. The consistent way of integrating the ry' yields L7 ~ 0{1) ^24;]. In summary, the Li in £4 are 0{Nc) at 
most, the Ci in Lq are 0{Nc) at most, and so on. 

Hence, if a diagram has Ngi vertices from Cd, they contribute, at most, with Nd{d — 2)/2 powers of Nc- Summing 
over all the d, the scaling of the LEG that contribute to that diagram is given, at most, by J2d^did — 2)/2 powers 

of Nc- Taking into account that the 1/F"- factors behave as 0{Nc^ ), we conclude that the non-factorization terms 

should be 0{Nc'^ '' ) = 0{Nc ) at most. But since we noted that non-factorization terms require 

1/ > 1, then the largest factorization breaking contribution is ©(A*""^), at most. Actually, this is the behavior of the 
non-factorization correction we explicitly calculated in Eas. ((29| and (|30)) . This ©(A""^) counting of the factorization 
breaking, which we have formally showed here in the low-energy representation, confirms what had been suggested 
previously in the literature ;2?^ . 

Finally, if we compare with the original QGD factorization hypothesis Eq.(IT|), we conclude that factorization of the 
four-quark condensate as the square of the two-quark condensate holds formally in the Nc —> 00 limit. This is of 
course only a formal statement, since we have just seen that in the low-energy calculation the factorization breaking 
terms diverge. 

VI. CONCLUSIONS 

In this work we have addressed the issue of the four-quark condensate factorization into the two-quark condensate 
squared, within the low-energy representation of those condensates provided by Ghiral Perturbation Theory. 

Our main result is the formal model-independent proof of the non- validity of the factorization or vacuum satu- 
ration hypothesis for the low-energy sector of QGD. A detailed calculation of the NNLO two-quark and four-quark 
condensates both for two and three flavors shows that, to that order, factorization is broken by terms which cannot 
be rendered finite with the usual renormalization procedure ensuring that the two-quark condensate is finite and 
scale-independent. This breaking of the factorization assumption at low-energies is then a mo del- independent result, 
since it relies only on the effective Lagrangian formalism, and is consistent with previous observations regarding the 
incompatibility of the factorization hypothesis with the QGD renormalization group evolution. In addition, the very 
same non-factorization term is obtained by using more conventional definitions of the quark condensate within the 
MS scheme in dimensional regularization. As a consistency check of our analysis, we have derived the light and 
strange susceptibilities from the calculated four-quark correlators, showing that they agree with a direct derivative 
with respect to the quark masses of the two-quark condensates. The explicit renormalized and scale-independent 
expressions for the GhPT NNLO susceptibilities are not given elsewhere, to our knowledge. Factorization would 
formally hold in the Nc ^ 00 limit, which we have been able to show to any order in the chiral expansion, the leading 
factorization breaking scaling as 0{1/Nc). 

We believe that these results can be useful for workers in the field, in particular concerning the OPE and sum-rule 
approach. A natural extension of this work is to consider finite temperature effects to see how they affect factorization 
and its connection with the chiral susceptibility, which in the thermal case plays a crucial role near chiral restoration 

M- 

Appendix A: Quark condensates to NNLO in ChPT and their renormalization 

In this section we will give our NNLO results for the two-quark condensates. As explained in text, the corresponding 
four-quark condensates cannot be obtained just by squaring these results, but one also has to add the non-factorizing 
contributions described in Eqs. (fT4| and (fT7|) . 

The free meson propagator in dimensional regularization is given by [8| : 

G,(0) = 2M^^-^X, (Al) 
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with 

r ri - ^1 

and Z? = 4 - e. 

The SU{3) £4 ChPT Lagrangian is weU known jS] and we do not reproduce it here. The relevant terms for the 
calculation of the condensates in the 0{p^) Lagrangian |l3| are only those dependent on the quark masses to leading 
order in the Goldstone boson fields. Here, we will follow, for simplicity, a different notation than in [lOl to denote the 
£q low-energy constants involved in the mass terms: 

^™„5a(3) ^ ^(^c.m^^C^m'm^+Csmml + C^ml), (A3) 

Recall that our Ci are linear combinations of the LEG considered in flO, '26\ whose precise form is not relevant here. 
Nevertheless, we still follow the convention in [2g] for the renormalization of the 0{p'^) and 0{p^) LEG in the MS 
scheme: 

k = m^-4[;[(m) + 7.a], 
h, = {c^lf-^[hl{^i) + s,A], 

L, - (c/i)^-4[L[(M)+r,A], 

H, = M^-Mi/[(/i) + rfA], 

a = ic^,r^-'^ [cKm) - f(^«) A^ - (ff + ff{p)) a] , (A4) 

where fi is the renormalization scale, A^^ ~ 16tt^{D ~ A), logc = — [log(47r) — 7 + 1] /2, 7 = — r'[l], 7^, Si, Ti, 

rf, 7j , Tl^'^ , 7I and F^ ' are numerical coefficients, whereas 7/^, 7^ are linear combinations of the L\{^). The 
above expression for the q shows that these constants have to absorb both two-loop divergences with £2 vertices and 
one-loop ones with one £4 and one £2 vertices. 

The renormalization of the Li in Eq. (jA4l) coincides with that in 8] up to 0(1) in the e expansion: 

i, = L[(Ai)+r,^^-*A + 0(e), (A5) 

and so on for the Hi, whereas the U^hi renormalization coincide with lll[ to that order. For the renormalization of 
the one-loop effective action, the 0(e) in Eqs.(jA4| and (jASP can be neglected. However, when two-loop diagrams 
are considered, as it is our case here for the quark condensates (e.g., diagram (d) in Figure [T]) products of the form 
LiG(O) yield finite contributions not vanishing in the e — > 0"*" limit. The 0{e) has to be kept also in the expansion of 
A in Eq. (|A2|) when expanding Gi(0) in Eq. (|Aip in Gi(0)^ contributions. 

As for the /.t scale dependence, the Li,li and the Ci,Ci are scale-independent so that the scale dependence of the 
L^{ij),ll{fi), CI {fj.) , cKfj.) is canceled with the exphcit /^-dependence appearing in Eq. (jA4|) . This allows to express 
all the logarithms of the masses referred to the scale /i, i.e., \og{Mf / ij,'^) so that the final result for the observables 
should be finite and scale-independent. 

We also recall that to the order we are calculating, the propagators are renormalizcd to NLO (tadpole corrections) 
and one has to include the wave-function and mass renormalization to that order. The renormalized masses are given 
in [8], while the explicit wave-function renormalization can be found for instance in |27|. We recall that we should 
include now up to 0{e) in those tadpole corrections, for the reasons just explained. 

With these renormalization conventions, we turn to the NNLO quark condensates. The F^ coefficients appearing 
in the calculation are [8|: 

73 = -1/2, ,5i = 2, 

F4 = 1/8, F5 = 3/8, F6 = 11/144, F7 - 0, Fg = 5/48, Ff ^ 5/24, (A6) 

Recall that in SU{3) Lg, Lg and H2 come explicitly from the £4 vertex contributions to the condensate and are 
therefore the only LEG appearing to NLO. The mass and wave function renormalization bring up also L4, L5 and L-j 
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to the final result. L7 only appears in the 77 mass renormalization. In the pure SU{2) case, only I3 and hi enter in 
the calculation. 

Once the above LEG renormalization is performed, we have checked that one can choose the Ci and Ci in Eq. (jA3|) . 
renormalized through Ea. (|A4p . so that the final result for the two-quark condensates is finite and scale independent. 
We obtain: 



7'^''' ^ 12, 



ri- 



7f = -48Z, 



3' 



\{sq) 



\{sq) 



1, r^"'^ = 32/27, rp' = 64/9, TX"" = 160/81 



(^?) _ 



32 



32 

r^ = y (162L^ + 31L^ - 324ig - 62i^) , 

32 

rf = —{96Ll + 35Ll~l92Ll + 2AL^T~62Ll), 

32 

ri = -{78Ll + 43Ll-6{26L', + 8L^r + 17Ll)), 



(A7) 



and all the linear terms 7} — F^- — for the above LEG. 

For convenience and following the same notation as [8|, we define: 



Mi 



— ^^^loK^, 



Vr = F' 






1 



In SU(2) the leading order pion mass is related with the physical one by 



/i2 



(A8) 



M^ = M^^ 1 + Ai. + 






(A9) 



and in SU(3), 



Ml =M? 



Ott 



1 , ^'f 

1 + A^TT - Y 



16Mo\^2ig - LD + ^ (2iS + 2i^ - L^ - Ld' 



F^ 



Ml =M? 



^K 



QK 



2jJ^ 

3 



8Ml 
F2 



l + ^ + ^^(2i6-i4) 



F2 



2Lg — 2L^ — Lr-^j 



M^ =m2. 



Mol 



1 + 2fj.K - -:;^^■r) 



8M0', 



QK 



-A^TT + ■7;^J■K 



;f^n 



p2 i2Ll-Ll) + ^{2M, 



Ml) {2LI LI) 



(AlO) 



The relation between the leading order pion decay constant and the physical one up to two-loops in given in [28| 
for SU(2) and in "2^ for SU(3). 

The final expressions for the two-quark condensates, finite and scale-independent, up to NNLO, that have been 
calculated previously in [1^] for SU{i), are given by: 



/- \SU(2) 
m)i.NLO 



\'i'ill,NNLO 



-2BqF^ i 1 



2AC 
F2 



ih[ + Q - 3m. 



{m) 



5(7(2) 
l,NLO 



2BqF^ 



-JK 



F2 



(MttJ^tt +4/3^^) 



^(-le^gV. + cD 



(All) 



(A12) 



12 



Wi/LNNLO 



-2BqF^ I 1 



^ [(i/J + 4Lg + 2LI) Ml^ + SigAfoV] " 3^. - 2^.^' - ^M, 



(99) 



S(7(3) 
l.NLO 



F2 

2BoF^ 






TMKMj? 



1 

F2 



3A'fo^/i^i^^ 



-A/f2 



OTrt^Trl^ri 



Q 4 1 

24 
+ -p;^/x. [i3Ll + 2LI - 6LI - ALl) M^^ + 2 (L^ - 2i^) Mq^] 

Ifi 
+ ;^A^^ [(^4 - 2iD ^^'. + 2 (3i^ + i^ - 6Lg - 2L^) M^V] 

-I2 Mr, [(-3i$ - 2i^ + 6LI - 48L; - 12LI) A^ + 2 (15L^ + 4L^ - 30Lg + 24L^) M^^] 



9i^ 

24Afo^ 

16M2 



F4 



-^. [(L^ + L^ - 2LI - 2LI) A/2^ + 2 (i^ - 2Lg) A^^V] 



[("3LI 



(A13) 



6L^-48L;-18L^)AC 



27J^4'^'7 

r 

4 

(acr ~ 2(72^ + q) Ml + 4 (C2^ - Cl) MlM^K + AC^,M^k 



2 {3LI - All - 6LI + 48L; + 24L^) Af ^^Af, 



2 

OK 



+ 8i3Ll + 2{Ll-3iLl + ^ + Ll)))M^^ 



1 

8F4 



(A14) 



{ss}nlo 
{ss)nnlo 



-^ [- {HI, - 4LI + 2LI) Ml + 2 (HI, + ALl + 2LS) Af^V] - ^^^K - ^M, 



Sof ^ \ l^^, 



-.IJ-KlJ-i^ 



1 
F2 



^Mlfi^i^,^ 



-BoF^ { 1 

(ss}a'lo 

+ f^^L^{L\~2Ll)Ml 

+ J2^'K [{LI ~ 2LI) Ml + 2 {2LI + L^ ~ 4ig - 2L^) A/qV] 

+ ^M. [(3il - 4L^ - 6LI + 48L? + 24L^) Ml + 8 (SL^ + 2(L^ - 3(L^ + L^ + L^))) A^^^] 



32Af, 



OK 



pi 

32 



i^K [(il - 2LI) Ml + (2LI + LI- iLl - 2LI) A/qV] 



(A15) 



■ 27i^4 ■^'J [(-3^4 + ^5 + 6i^ - 48L; - 18L^) Af ^^ + 2 (3L!; - 4L^ - 6L^ + 48^^^ + 24Lg) A^ol^^oV 
+ 8 {iL\ + 2(L?i - 3(Lg + L; + L^))) AfoV] 



+ 



1 

4F4 



(4^ - 2Cl + ici) Ml + 4 {Cl - 3Cl) MlM^.j, + UC'IM, 



4 

OAT 



(A16) 



where the Gell-Mann-Okubo relation 3Mq = 4AfQ^ — Ml for the SU{3) leadmg order masses has been used and 
the renormalized L[, l^ and c[, Cl constants depend on the scale 11 as explained above. 

Appendix B: Four-quark condensates in the usual AiS definition 



Here we consider the definition in Eq.(|31) of the four-quark condensate in Euclidean space. Let us restrict to SU{2) 
since it will become clear that the argument can be straightforwardly extended to the S'C/(3) case. The four-quark 
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correlator to NNLO is given in Eq. p4)) . so that its Euclidean Fourier transform to this order is (see our Euclidean 
space-time conventions in section llll| : 

n(Q2) ^ (2^)^(gg)25(^)(g) + 2Bl [A{h + h,) + 3J^{Q^)] (Bl) 

with Q2 ^ ^^^^ Q2 and 



(27r) 



MQ^) = / 7KzmG-(K)GAK - Q), (B2) 



which is nothing but the one-loop integral appearing in pion-pion scattering, dimensionally regularized in [ll|. Its 
divergent part is contained in J^(0) = -2M^-'^X - l/(167r2) with A defined in Eq.^M^ while J(Q2) = J^{Q'^) - J^(0) 
is finite. Note also that J-r^iCp) defined in Euclidean space is real. The imaginary part in J.^ giving the usual unitarity 
cut in scattering amplitudes arises when the analytical continuation of Q^ to Minkowski space-time is performed, but 
here we should keep the Euclidean version, since we are following the prescription in Eq. ^ to perform the additional 
momentum integral. 

Before proceeding to the calculation of the four-quark condensate, let us note that the divergent part of the J^ in 
Eg. dBip cancels exactly with the LEG contribution since l-^ + hi — l^ifJ.) + hKfj.) + (3/2)^^~'*A (see Eas. (IA4l) and 
(jA6P ). Thus, n((3^) is finite and scale-independent before integration in Q. This is actually a welcomed check, since 
the scalar susceptibility given in Eq.dm can be written also as xi = n(0) with n(Q^) = U{Q^) - {2Tr)° {qq)^ S'^^'^ (Q) 
and should be finite and scale-independent. 

However, we will immediately see that the additional integration in Q in Eq.Q generates an extra divergence which 
cannot be removed and in the end gives the same divergent factorization-breaking result as the definition in Eq.Q. 
For that purpose, let us follow the standard dimensional regularization procedure [22| and write: 

•^-(Q') = 7rTD72 / '^^ / dW'^''/^ exp {-A [M^ + Q^xil ~ x)] } (B3) 

(,47rj ' Jo Jo 

which is valid within the domain Re [D] < 4. Now, before performing the x and A integrals above, we integrate over 
Q so that: 



MQ')^7TTn\ dx[xil~x)]-''^' 



{2-k)D "^- ' (4^)^ Uo ^ ^ " J Uo (4vr 



= G^(0) 



(B4) 

where the one-dimensional integrals are solved for Re \D\ < 2 and we have used standard properties of the Gamma 
function. Since the result is analytic in D, it can be extended to Z? = 4 — e with e -^ 0+. Therefore, integrating in 
Ea. (IBl[) over Q according to Eq.Q and taking into account that J d^Q/{2'K)^ — (5'-^'(0) = 0, gives exactly the same 
divergent factorization-breaking result for the four-quark condensate as the one using the prescription of Eq. ([2]) . 

Another way to arrive to the same conclusion is to perform the change of variables Q -^ Q + K in the double 
-D-integral J d^Q J d^K in the region of D where it converges, which in this case is Re [D] < 2, which follows by 
direct power counting in Q and K of the propagators in Eq. (jB2[) in the large Q^ and K^ Euclidean region. 

It is clear that the same equivalence between the two definitions holds in the SU{3) case simply by considering Jk 
and Jjj apart from Jn-, since the results of the correlators in Egs. p^ and ([20]) do not mix different meson species. 
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